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SUMMARY 


This report documents the results of work accomplished by the authors on 
a task undertaken to suggest an appropriate approach to the solution of the optimum 
two-impulse transfer problem between orbits of specified inclination. The task in- 
cluded a literature search to identify the current state of the art and a definition of 
the suggested approach for the specific application of a lunar orbit trim. Some num- 
erical data available at the start of the study prompted Mr. Pines to undertake some 
original research on the subject which was outside the scope of this task. Although 
the research study is not reported here, the applications of his results to the pro- 
blem of interest is included as Appendix A of this report. The formulation for a 
computer program developed under this task following a more conventional approach 
is included as Appendix B, 
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I. The Problem 

The first phase of this assignment was to conduct a literature search to 
discover what numerical and/or analytic techniques might be available for the 
solution of the following impulsive orbit transfer problem: 

1. All orbital parameters are given for an initial elliptic orbit. The 
only assumption made on these parameters is that the eccentricity has an upper 
bound of .8. 

2. The final orbit is to be circular, have specified inclination and 
radius specified to within ± 200 km. 

The application is to a circular orbit about the moon at critical inclination. 
The initial orbit will be achieved by a retro -burn from the hyperbolic approach 
trajectory, and will be tracked long enough to establish its elements. One would 
hope that this initial orbit would be close to the desired final circular orbit, but 
techniques are desired for achieving the optimal insertion into the circular orbit 
whether or not this hope is realized. 

The second phase of the assignment was to pursue analytic investigations 
suggested by the literature search. 

n. Results of the Literature Search 

The starting point for the search was a comprehensive survey by Gobetz 
and Doll^ published in May 1969, which covers work on impulsive transfers 
carried out before mid-1 968. An extensive bibliography (over 300 references) 
is included, and an outline of general results obtained is presented. 

An effort was made to review all of Gobetz' s references which appeared 
at all pertinent to the problem at hand. It was not possible to obtain copies of 
all such references; however, enough were reviewed to be reasonably sure 
that no general analytic conclusions are available in the literature for this pro- 
blem - as, for example, the conclusion that if the orientation of one of the orbits 
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is left entirely unspecified, optimality of the transfer requires that this orbit be 
chosen coplanar and coaxial with the other orbit. 

The papers actually examined from Gobetz's list of references are those 
numbered as follows: 

9 , 17 , 19 , 20 , 35 , 37 , 54 , 55 , 

56 , 61 , 104 , 105 , 106 , 113 , 123 , 

124 , 125 , 192 , 193 , 195 , 213 , 

214 , 216 , 220 , 221 , 241 , 251 , 

276 , 277 . 

In addition, recent issues of Astronautica Acta, the AIAA Journal and the Journal 
of the Astronautieal Sciences were scanned, and several additional references were 

/Q \ _ /Q\ 

examined ' ' * . Prior familiarity with the work of the groups at Huntsville 

(Hoelker et. al. ), G. E. (Altman et. al. ), North American (Bender, McCue et. al.), 
Stanford (Breakwell and his students) and the work of F. T. Sun made it unnecessary 
to look in detail at all of the papers by these authors listed by Gobetz. 

It should be noted that my examination of these papers was primarily 
directed at looking at the results claimed by the authors in order to determine 
applicability of their work to the particular problem outlined above. No effort 
was made to follow all the analyses through in detail. Of all these papers, very 
few dealt with non-coplanar transfers and still fewer with transfers of more than 
two impulses. No useful analytic results w r ere found. H. W. Small T s program 
for optimal N-impulse transfers appears to have possibilities for our problem. 

This program is a culmination of the efforts of Breakwell and several of his 
students over some years. It is based on the indirect method, requires some 
preliminary numerical work, and requires guesses on the initial values of the 
adjoint variables. All trajectories generated are optimal and the iteration in- 
volved is to match the terminal conditions. Small has used the program for an 
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exhaustive study of coplanar transfers. Our problem is not, however, the pro- 
blem of surveying all possible transfers, but rather of optimally matching certain 
terminal conditions. S. Pines has criticized procedures such as this on the ground 
that the numerical match obtained on the terminal conditions will not be adequate 
for practical implementation. Small has made considerable additional material 
(in the form of a very rough draft of his thesis) available. Some effort would be 
involved in understanding the procedure, if one wished to use his program on this 
problem. Its principal attraction appears to be that it would yield the global N- 
impulse optimum. 

m. Theoretical Work 

Construction of a direct theory for a two-impulse solution to the problem 
was started, but the time limitations did not permit its completion. The para- 
meters tentatively selected were: 

1. Argument of latitude of the departure point on the initial orbit. 

2. A parameter Cl related to the radius of the circular orbit by 

r = a + p cos ct a given 

p - 200 km 

3. Longitude of the ascending node of the circular orbit. 

4. Argument of latitude of the arrival point on the circular orbit. 

These four parameters define the initial and final position vectors on the 
transfer orbit, and one or two transfer- orbit parameters are required as follows: 

5. Any inplane parameter of the transfer orbit: semimajor axis, 
eccentricity, argument of latitude of one end or the other, magnitude of the angu- 
lar momentum, or transit time. 

6. If the initial and final position vectors of the transfer orbit are 
collinear, an out-of -plane parameter: either inclination or longitude of the 
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ascending node. 

The formulation of the problem is such that any choice of parameters which 
define the geometry of the transfer is easy to implement. As noted earlier the 
theory is not complete; further detail does not seem warranted for reasons which 
follow immediately. 

About midway through the study, the conjecture (partly based on numerical 
results of C. Uphoff) was put forward by S. Pines, that the optimal two-impulse so- 
lution to the open-time problem is a 180° transfer chosen such that the changes in 
angular momentum resulting from the two impulses are collinear. The proof of 
this will be published independently; however, the application of this solution to the 
problem at hand is outlined in Appendix A. 

To permit the study of fixed transfer-time transfers, a two-impulse tra- 
jectory optimization program based on the indirect technique was constructed. 

A general iterator was incorporated in the program which also permits optimi- 
zation by direct methods. Although the formulation for this program is suffi- 
ciently general to include open- time problems, it will not yield the 180° transfer 
(which appears to be the globally optimal solution) because such a transfer re- 
presents a singular solution for the formulation employed. The formulation is 
described in Appendix B. 

IV. Conclusions 

1. The general two-impulse trajectory optimization program is 
available for fixed transfer- time problems. 

2. Pines 1 two-impulse 180° analytic solution is available for open- 
time problems. 

3. For multi-impulse solutions, SmalPs program is a possibility. 

As noted above, considerable preliminary study is required. The detailed 


5 



material supplied by Small did not arrive in time to carry out such a study, nor 
were computer facilities available for familiarization with the program. 

4. Because of the last comment, no assessment of the relative merits 
of Small’s program compared to the first two items will be attempted. 
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APPENDIX A 


OPTIMAL TWO-IMPULSE TRIM FOR RAE-B 


INTRODUCTION 


The equations for the optimal two-impulse trim for transferring from an 

elliptic conic with given a , e , i , O , W , and t to a required circular orbit of 

P 

given inclination are described herein. The final circular radius is designated by 
r . A test is included for the condition where a nonoptimal three- or four-impulse 
trim may prove superior for conditions of large changes in the angular momentum 
vectors. The results are given without derivation. 

NOTATION 

i initial inclination 

o 

i^ required final inclination 

H q initial angular momentum vector 

A 

k unit vector in the direction of the equatorial pole of the planet 

h magnitude of H 

o o 

h^ magnitude of the final circular angular momentum vector 
unit vector in the direction of the initial ascending node 
r* 180 ° transfer radius on original conic 
r” radius opposite r+ on original conic 

h^ magnitude of the transfer angular momentum vector 
eccentricity on initial orbit 
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DESCRIPTION OF THE SOLUTION 


The optimal two-impulse solution is characterized by a 180° transfer. 
The direction of the two-impulses are anti -parallel. The initial, transfer, and 
final angular momentum vectors all lie in a plane and the two vector changes, In 
angular momentum are collinear. 


EQUATIONS DEFINING THE SOLUTION 


Let the minimum possible change in inclination be Ai^, then 


Ai = i. - i 
M f o 


The initial unit ascending -node vector is given by 


n = 


kxH 

c 

sin i 


where 


H * 


H /h 
o o 


The minimum possible change in angular momentum, corresponding to 

^ A 

Aij^., is given by rotation of H q about the initial ascending node through 

the angle Ai,,. 

M 


fM 


H_ = h,£^cosAi„ + 


sin Ai„, cos i 


sin Ai, 


M 


sin i Jo sin i J 


where 



?r> 


Introducing the single-degree-of-freedom variable, X, as the angle be- 
tween the initial and final ascending nodes, we have 

A - ~ - 

0. = cos X £} + sin X kxO (4a) 

to ' ’ 


The final angular momentum vector is given by a rotation of about 

through X. 


H. = cos X H-__+ (1 - cos X) cos L k + sin X kxH„, 
f fM ' 7 f fM 


(4b) 


The corresponding increment in inclination from the initial angular momentum 
vector to the final angular momentum vector is given by 


, H xHJ 
A * . “I O f' 

AL = sin — r— r — 

f h h, 

o f 


(5) 


The 180 transfer takes place along the unit position vector on the original 
conic and is given by 


R = 


H x H, 
o f 


h h. stn A l- 
of f 


( 6 ) 


The central angle between the original ascending node and R is given by 


-1 A 

cl - cos R) 


(7) 


The conic radial distance in the positive direction of R is given by 


+ _ o 

*1 fj, £ 1 + e^ cos ( a - ttT)J 


( 8 ) 
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The radial distance along the negative direction of R is 

, 2 


r = 


^[1- e Q cos(a-tO o )] 


( 9 ) 


The first impulse takes place along the larger of the two radii. If r*>r^ » then 
the radius vector for the first impulse Is 


R i= r i +R 


(10a) 


and the second impulse takes place at the radius vector 


R 2 “ r r 


(10b) 


If r "" > r^ + , then 


R i “ - r i' R 


R 2 = r R 


(U) 


Figure 1 is convenient to describe the various angles and momenta mag- 
nitudes required for the solution. 
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Let the vector change in angular momentum be AH 


AH = H. - H 
f o 


and its magnitude be 


6h 


•/ 


2 2 
h . - 2 h, h cos A i _ + h 
f f o f o 


From the law of sines, we have, for the angle between H and AH, 

o 


f 

sin a - — sin A i, 
h f 

o 


and 


cos a = 


h - h, cos Ai, 
o f f 

6h 


The transfer angular momentum has a magnitude given by 


_ / 2 ^< rr i + ) 

“t 4 - + 

r + r. 




or 


h = 
T 


2/i( r < 


if r> 


r + 


■i 


In either case, we have, for the angle between H,^ and AH, 


h 

o 


sin /3 = sin a — 


( 12 ) 


(13) 


<14a) 


(14b) 


(15a) 


(15b) 


(16a) 


A -S 



(16b) 


cos |8 = sign ( hy, - h Q cos a) J X - si 


sin^ 0 


The angle between H and is 

o T 


Al t - TT-(a+l 3) 


The transfer orbit angular momentum vector is given by a rotation of about 
the unit vector R^ /r^, where r^ is the larger of r^ + and r^“ , through the 
angle Ai T , 

H T = {” cos (“ + 2> V Sln < r a+ ^ ) BjXhJ 


(17) 


The magnitude of the change in angular momentum after the first impulse is 
given by 


6h = |H_-H | = — r^- [sin (a+0)] 
1 T o l sina 


(18) 


Let the vector dot product of the position and velocity prior to the first 
impulse be 


d = R • R 
1 1 1 


(19) 


Then, the vector dot product after the first impulse Is given by 


<- 


6h. 


1 - 


(Sh-e^) 


( 20 ) 
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The velocity vector after the first impulse is given by . 



h t xR i d l +R l 
— — 


The velocity prior to the second impulse Is given by 



H vR A 
T 2 

-2 

r 


< H 2 

-2 

r 


The final circular velocity vector after the second impulse is given by 



h t *r 2 




We may now compute the sum of the two impulses and carry out the one- 
dimensional search for the X which yields the minimum sum of the two impulses. 

6v i = l R i- R i'l 

6V 2 = 1 R 2- R 2 'I 

This is the optimal two-impulse trim unless the condition in the next section is 
violated. 


TEST FOR MULTI-IMPULSE NONOPTIMAL TRIMS 


More than two trims will be better than the two-trim program outlined 
herein if 



( 26 ) 


H(/^ -Olv^ -O h 


For such an event, there is no unique minimum and much better three- or four- 
impulse trims can be found for longer transfer times. 



\ 
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APPENDIX B 


GENERAL FORMULATION FOR OPTIMAL 2-IMPULSE 
TRANSFER BETWEEN INCLINED, ELLIPTICAL ORBITS 


Consider a spacecraft in an initial orbit with elements a , e , i , O , w , 

o o o o o 

and t which denote the semimajor axis, eccentricity, inclination, longitude of 
po 

ascending node, argument of periapse and time of periapse passage, respectively. 

Suppose it is desired to transfer to a final orbit of specified semimajor axis a^, 

eccentricity e^ and inclination i^ using two impulses, the magnitudes of which 

we denote v and v., so as to minimize the sum It = v + v*. The motion is des- 
o f of 

cribed by the vector equation 

™ 3 

R - - ^ R/r 


which represents free flight in an inverse -square force field where R is the 
position of the spacecraft, r= | R | , and \i is the gravitational constant. Clearly, 
the position at the time of the first impulse R^ and the velocity before adding 
the impulse R q are functions only of the time at which the impulse is applied, 
t.e. , 


and, thus 


R = R (t ) 
o o o 


R = R (t ) 
o o o 

R + = R ~(t ) + V 
o o ' o o 


where the plus sign denotes the limit evaluated just after the impulse and V q is 
the impulse velocity vector. With a^, and i specified, the position and 
velocity after insertion in the final orbit may be considered to be functions only 
of the final longitude of ascending node and the argument of latitude jS^, i.e.. 


R f = R f (O f . /3 f ) 
R f + = R f + (O f , /3 f ) 
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and therefore 


R f " = R f + (Q f , j8 f )- V £ 

where is the final velocity impulse. 

Proceeding formally with the application of the Maximum Principle, one 
may write the adjoint equations for the problem which define the behavior of the 
primer vector A, 

A ~ (R * A)R - A 

5 o 

r r 

This equation is known to possess the analytic solution 

where <§(t , t) is the state transition matrix for the two-body transfer from 
o 

time t to time t. The transversality conditions for the problem are obtained 
o 

from the Maximum Principle by applying the general equation 

dff + [A • dR - A * dR - hydt] o = 0 

M . * • 

where h = A - R - A • R is the variational Hamiltonian which is a constant of 

the motion. Note that the differentials dR are evaluated at the terminals of the 

• +* — 

transfer trajectory; hence, dR Q and dR f must be used. Forming the appro 
priate differential, one obtains 

• 

dR = R ’ dt 
ooo 

dR + = R dt + dV 
o o o o 
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dR = (k x R ) dfJ { + (h x R ) d/3 

dR f " = (it x R f ) dO f + (ii f x R f + ) d/3 f - dV 


where k is a unit vector in the direction of the North pole of the central body and 

V is a unit vector in the direction of the angular momentum of the final orbit. Let 
f 

each of the velocity impulses be defined in terms of a magnitude and two angles 
which uniquely define angular orientation. We arbitrarily select one angle, a, to 
be a rotation about the (initial or final) unit angular momentum vector h and the 
other, y , to be a rotation about the unit vector j defined by j = (V x h)/|v x h|. 
Then 

V 

d V = — dv + (h x V ) da + (j x V ) dy 
ov o oooooo 
o 

v f 

dv. = dv, + (h x V ) da + (j x V ) d y t 

1 Vj f f f f f f f 

Finally, writing 


3 TT 5 77 

dTT — T"“” dv + — dv x = dv + dv r , 
ov odv. f o f 

o f 

substituting into the general transversality condition, collecting coefficients of 
the remaining independent differentials and setting to zero, one obtains the follow- 
ing transversality conditions 


(a) k • (F + V f x A f ) = 0 

(b) y (F+V f x^) = 0 


(c) 1 - (A f • V f )/v f = 0 

(d) h f '(V f x^) = 0 

(e) j f • (V f x A f ) = 0 
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(f> 

1 - (A • V )/v = 0 
0 0 0 

(g) 

- h • (V x A ) — 0 
0 0 0 

(h) 

I 

1-1.1 

0 

■ 

> 

0 

II 

0 

(i) 

0 

II 

> 

rC 2 

0 ) 

- A • R + A * R 
000 


where F is a vector constant of the motion defined as follows 


F=RxA-RxA 

The equations (d), (e), (g) and (h) along with the definition of j are sufficient to 

conclude that the optimal directions of V and V, are collinear with A and 

of o 

Aj, respectively. Equations (c) and (f) indicate that the velocity increments 
must also have the same sense as the associated primer and that 

|A o l=|A f |=l. 

Due to the collinearity of the velocity impulse and primer at the final time, equa- 
tions (a) and (b) reduce to 

(a') k * F = 0 

(b 1 ) h f • F = 0 

which means that F will be aligned with the line of nodes of the final orbit. 
Equations (i) and (j) arise if final time and initial time, respectively, are 
left completely open. If, instead the flight time is fixed, then the two equations 
(i) and (j) are replaced in favor of the single equation represented by their 
sum, i.e.. 
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(j r ) 





= 0 


The above formulation may then be implemented in the following manner. 


(1) Estimate the initial time, final time and the longitude of ascending 
node of the final orbit and the argument of latitude at insertion into 
the final orbit. Along with the specified orbital parameters, these 
are sufficient to define the state and time in the initial orbit at de- 
parture and in the final orbit at arrival. 

(2) Solve Lambert's problem for the transfer trajectory between the two 
terminals, Simultaneously, evaluate the state transition matrix for 
the transfer. 


(3) 


Evaluate the velocity impulses by differencing the velocities at the 
terminals of the transfer orbit and the corresponding velocities in 
the initial and final orbits, i. e, , 


V 

o 




• + • — 



(4) Define the initial and final primer vectors 


A = V /v 
0 0 0 


VY v f 


(5) Partition the state transition matrix into four 3x3 elements 



and solve for A and A. as follows 
o f 

A = B^iA. - AA ) 
o ' f o . 

A = CA +dA 
f o o 


and evaluate the constants of the motion F and h^ . 

(6) Solve for the left-hand side of equations (a 1 ), (b r ) and (i), (j) or 
(j ! ) and At. 
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(7) Initiate an iteration of steps (1) - (6) to solve the boundary value 

problem in which t , t^, and are the independent parameters 

and the parameters evaluated in step (6) are the dependent para- 
meters. 

A computer program which successfully implements the above steps has 
been developed and delivered to the Principal Investigator. Limited exercising 
of the program prior to delivery indicated that the approach exhibited strong con- 
vergence properties. A major limitation of the program is that it cannot handle 
a 180° transfer trajectory which appears to be the globally optimal solution for 
an open-time transfer. The open-time solution which the program yields is a lo- 
cally optimum solution which, for the specific cases investigated, yielded a velo- 
city impulse sum about 30-40 percent greater than the 180° transfer solution. 

Care should also be exercised in solving for fixed-time solutions with the program 
because the existence of multiple local solutions to such problems are common. 

The particular local solution which the program converges to is totally dependent 
upon the initial guess of the independent parameters. Therefore, if the existence 
of multiple solutions is suspected, one should input estimates of the independent 
parameters in the neighborhoods of each suspected solution to find the global so- 
lution. The alternative is to map solutions throughout the admissible ranges of 
the independent parameters. 
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